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Abstract. 

Yield Optimized Interpolated Superoscillations (YOIS) have been recently 
introduced as a means for possibly making the use of the phenomenon of 
superoscillation practical. In this paper we study how good is a superoscillation that 
is not optimal. Namely, by how much is the yield decreased when the signal departs 
from the optimal one. We consider two situations. One is the case where the signal 
strictly obeys the interpolation requirement and the other is when that requirement is 
relaxed. In the latter case the yield can be increased at the expense of deterioration 
of signal quality. An important conclusion is that optimizing superoscillations may be 
challenging in terms of the precision needed, however, storing and using them is not 
at all that sensitive. This is of great importance in any physical system where noise 
and error are inevitable. 
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1. Introduction 

Superoscillatory signals are band limited functions that oscillate over some regions with 
a frequency larger than that of its maximal Fourier component. A number of examples 
have been given in the past for such functions with very interesting applications in 
quantum mechanics [Il[2l[3llll[5l[6l|7l|8], signal processing [9], [101 IHl [I2l ESI EH ES] and 
in optics, where superoscillations are intimately related to super-resolution [IHIETIEE] 
with some recent exciting experimental achievements [T^ |20] . For a recent review see 
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A key issue in the field is quantifying the superoscillations using some local 
frequency, for example by taking the local rate of change of the phase of the signal 
ra among other options. Another question of much interest is how superoscillations 
interact with physical matter during propagation [23l | 23 ] , with some indication (both 
theoretically es and experimentally m) that nonlinearities may actually play an 
important role in utilizing them. Furthermore, ref. |2^ provides a prediction of a 
quantum effect directly related to superoscillations, causing momentum ’superkicks’ 
imparted to a small particle located near optical vortices. 

However, although the phenomenon of superoscillation suggests many practical 
applications, its uses are presently rather limited, because the superoscillatory signal 
exist in limited intervals and that the amplitude of the superoscillations in those regions 
is extremely small compared to typical values of the amplitude in non-superoscillating 
regions [a [la US]. This often renders algorithms for generating them numerically 
unstable [28]. One of the methods for obtaining superoscillations that lends itself easily 
to improving the situation described above, is that of forcing an interpolation within 
the prescribed interval that will necessarily lead to the required superoscillation [131129] . 
Actually, fitting functions to prescribed forms inducing superoscillations by interpolation 
has been applied to optical superresolution in [30ll3T], and the resulting matrix inversion 
problems were explored in [32] . 

We start our discussion by explaining the concept of interpolated superoscillation 
[131129] . Consider the signal 


fit) 
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which is periodic with a period of 27r. Choose a sub-interval [0, a] with a < tt and impose 
on the function f{t) M constraints in the interval, such that f{tj) = fij for 0 < tj < a 
and j = 0,..., M — 1 The constraints result in a set of M linear equations in -|- 1 
unknowns of the form, 
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where 


1 j cos{mtj) m 7 ^ 0 

7 ^ m = 0 • 

Generically this set of equations: has no solution for M > 
M = N + 1 and a whole space of solutions for M < N + 1. 
interest [131129] is 


( 3 ) 

iV -|- 1, has one solution for 
A particular choice of great 


a 

^ M-1^ 


and /ij = (—!)■’. 


( 4 ) 


Provided M < -|- 1, this choice constrains the function to oscillate within the interval 

[—a, a] between the values ±1 with a frequency 
7r(M- 1) 


u = 
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( 5 ) 
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It is thus clear that the frequency of oscillation within the interval [—a, a] can be 
increased indehnitely regardless of the largest frequency cos{Nt) included in the signal 
dO. As is well known, this comes, at a cost. The energy in the superoscillatory region 
is extremely small relative to the total energy. Therefore optimization is required if the 
use of the phenomenon of superoscillation is to become practical. 

A natural step is to optimize the superoscillating function for hxed a and M < A^+1. 
However, one has to decide hrst in what sense to optimize it. One approach [T3l [29] 

OD 

is to consider the energy of the signal, E = f p{t)dt, use the fact that / is band 

— OD 

limited and minimize the energy under the interpolation constraints. In the periodic 

TT 

function case, the equivalent would be to minimize, E = f p{t)dt under the choice of 

— TT 

the Fourier coefficients given by equation ([1]) and the constraints (jlD . A step forward 
would be to maximize the superoscillation yield, defined as 

/ f{t)dt 

Y{M,a) = ^ -, (6) 

/ P{f)dt 

— OO 

rather than the total energy. As will become evident in the following the two methods 
of optimization are very close to each other, provided the optimal yield is indeed very 
small, which in many situations is the case. It is in place to say that defining the yield in 
terms of the energy of the signal is one option, which relies on the function being square 
integrable, foe example. However, square integrability is not always relevant because it 
can be related to features of the signal far away from the superoscillatory region, such 
as a slow decay of the tail, see |33] for such an example. Therefore, more generally, 
one can devise other definitions of the superoscillatory yield which are more localized 
around the superoscillatory region. 

Plugging a signal described by its Fourier components, as in Eq. CDi, into the 
definition of the yield, given by Eq. (I6|), we can obtain the following expression for the 
yield 
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where the entries of the matrix A are given by 
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To implement the optimization of the yield [22], the original degrees of freedom, are 
rotated, going from the original set of the A’s to a new set of degrees of freedom B, as to 
break the linear iV+1 dimensional space, into a direct sum of two subspaces. The hrst M 
dimensional subspace is the space, S'm, where the M coefficients {Bn} are independently 
constrained and attain prescribed values determined by the constraints. Thus, from all 
that subspace, we choose a specihc vector B^. The remaining [N +1 — M)-dimensional 
subspace, S'jv+i_m, is totally unconstrained. 

The set of interpolating signals is dehned thus by 

B = Bat+i-m © , (9) 

where Bn+i_m G 

In Ref. [29] we have optimized the yield under the interpolation constraints. In this 
paper we explore the effect of noise on the superoscillatory signal. Such noise is abundant 
in many physical systems and must be studies in order to understand the applicability of 
superoscillatory signals. This question can be reformulated in lay terms: by how much 
is a typical superoscillating signal worse than the optimal? In section 2 we discuss the 
effect of noise on the distribution of the yield in the subspace of interpolating signals. We 
consider two common models of noise, namely a bimodal and a Gaussian distribution, 
and show that small deviations do not affect the yield, while larger deviations result 
in a yield which is orders of magnitude smaller, with a narrow transition zone between 
these two regimes. In section 3, we consider deviations from the Fourier components of 
the optimized superoscillatory signal in real space. Such deviations affect the quality of 
the interpolation, and thus may damage the superoscillatory nature of the signal. We 
therefore dehne and study in section 3.1 the interpolation sensitivity of the signal, which 
quantihes the extent to which the signal deviates from the interpolation constraints. 
We then explore in section 3.2 the impact of deviations from the Fourier components 
of the optimized signal on the superoscillatory yield. An interesting effect is that a 
moderate degree of noise can actually increase the yield, at the expense of modifying 
the shape, which could be of interest in applications. The last question we address in 
section 3.3 is the effect of these deviations on the number of oscillations that appear 
in the superoscillatory region. It turns out that large enough deviations can sometimes 
lower the number of oscillations rendering these signals not superoscillatory anymore. 
In section 4, we summarize the paper and draw conclusions. 

2. Noise in the subspace of interpolating signals 

In order to progress, we consider the functional dependence of the yield on the 
constrained and on the unconstrained R’s. This composition is useful because any 
choice of the unconstrained R’s yields a signal which obeys the constraints. 

The numerator I on the right hand side of equation ([7]), can be expressed as follows 

N N 

1= Y, ^mnAmAn = ^ 

m.n=0 m,n=0 
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where = RAR ^ , R being the rotation that takes the coefficient vector A into 
B. Let us describe the matrix A^^^ by the following block form 
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(11) 


where is the transpose of F. The yield can be now expressed in terms of the 
unconstrained B's as 
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To simplify the notation we rewrite 
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^mnBmddn T 'f'', T dDi 
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where the symmetric positive dehnite matrix A, the coefficients and the constants 
Di and D 2 depend on the interpolation constraints (and also on the specihc frame of 
reference used in the Sm subspace). 

We would like to take a random set of R’s and get an idea of the corresponding 
random yield to which they give rise. Such randomness serves two purposes. First, in 
many physical systems noise is prevalent and it is of great importance to understand its 
impact on the superoscillatory signal and in particular on its yield. Second, it provides 
some insight into the orders of magnitude of R’s which are anticipated in the context 
of interpolating signals. 

It can be shown that 0 ;^, Di and D 2 are typically much smaller than the quadratic 
terms in Eq. flT^ . This means that for a random choice of the R’s the yield would be 
typically given by 


Y ~ 


N-M 

E 

m,n=0 


^mnBmBn 
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y 52 

Z—/ m 
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(14) 


Thus, it is clear that typically, the value of Y is distributed between the largest and 
smallest eigenvalue of A, which spans many decades (the condition number of A, 
dehned as the ratio between the largest and smallest eigenvalue is very large, and so the 
problem is typically ill-conditioned). To give just two examples, the largest and smallest 
eigenvalues of A (^max, Ein) for the cases of A = 6, M = 4 and N = 9, M = 6 are given 
by (0.000179,3.1210“^'^) and (4.6510“®, 3.6410“^^) respectively. It is interesting to note 
that the corresponding optimal values of the yield are Y = 0.0138 (for N = 6, M = 4) 
and Y = 0.0004 (for N = 9, M = 6). This implies that a random choice of a signal that 
interpolates the constrained points gives a yield which is orders of magnitude less 
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than the optimal one, which emphasizes the necessity to optimize the snperoscillatory 
signals. 

At this point we can say more abont the statistical properties of the yield of 
interpolating signals. We will assnme hrst that the i?’s are independent and drawn 
from a distribntion 

P(B) = A (f) . (18) 

where g{x) is a normalized and symmetric distribntion, with variance eqnal to 1. 
Therefore, a measnres the width of P{B). Snch noise describe possible sitnations where 
the optimized set of i?’s are corrnpted to a certain extent, or that their accnracy has 
been degraded. Independence of the 5’s simply assnmes that the a random bias has 
been introdnced directly in the snbspace of interpolating signals. Namely, snch noise 
will not modify at all the fact that the signal generated by this set always obeys the 
snperoscillatory constraints. 

We will consider averages as a fnnction of a and then distribntions for hxed a. We 
start with some obvious results: it is clear that 

lim Y = — . 

a —^0 

This gives, combined with equation (1TT|) . a lower bound on the optimal yield, namely 

Yopt ^ max ^ ^ , l^max^ • 

Consider next the very simple average of the total energy per period (the denominator 
of Y eq. (E])) 

(E) = (iV + 1 - M)a^ + D2 , (18) 



which turns out to be independent of the specihc distribution. The average yield 
is however more difficult and does depend on the distribution. Let us therefore 
discuss a two illustrative examples, namely a bi-modal distribution, which is a discrete 
distribution, and a Gaussian distribution, which is continuous. These two distributions 
represent a range of possible symmetric distributions with a finite variance, which are 
often used to model noise. 

Consider first the bimodal distribution. 


P{B) 




(19) 


where 5{x) is the Dirac delta function. The average yield in this case is simply 


(>^) = 


tr A ] + Di 


( 20 ) 


(A + 1 - M)a2 + D 2 ■ 

This is nothing but a Sigmoid function, attaining the value D 1 /D 2 when a —)■ 0 (in 
agreement with Eq. (1T6|1 . and trA/{N + 1 — M) for a —)■ cx). Actually trA is very small 
for snperoscillatory signals [29], and thus (Y) is monotonically decreasing, interpolating 
between D 1 /D 2 and trA/{N + 1 — M) 1. 
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Figure 1. The average yield (Y) as a function of cr for TV = 6, M = 4. The upper 
curve is for the Gaussian distribution (given by eq. (1221) 1 and the lower one for the 
bimodal distribution (given by eq. (I^ l. 


Consider next a Gaussian distribution, 

Although considerably more complicated, the average yield can also be obtained 
analytically here. 




D, 

2a2 
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20-^ E N+l-M 
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where Ey (z) = f dx , is the exponential integral function [35]. In order to obtain 

1 

more insight, we can study the various limits o this expression. In fact, Ey{z) ~ e~^jz 
for 2 ) —)■ oo and therefore, it is easy to verify that here too (F) —)■ D 1 /D 2 in the limit 
(j —)■ 0 (consistent with Eq. ffT6|) b In the other limit, namely a —)■ cxo, we can use 
Ey{z) ~ l/(z/ — 1) for z —)■ 0 (and n > 1). Since N > M + 1, this gives rise to 
(Y) —>■ tr A/(A^ + 1 — M). Also, by differentiating Eq. fl2^ with respect to a it is easy 
to show that it is a monotonically decreasing function, since dEy{z)/dz = —Ey_i{z) 
which is negative definite for all real values of . 

It is clear from equations o and ffTbjl that regardless of the distribution the 
average has to interpolate between D 1 /D 2 at small a and trA/{N + 1 — M) (which 
is known to be a very small number [29]), and are monotonically decreasing. The two 
averages, (12(1 and fl2^ are presented in Figs. [Hand[2] for the cases N = 6, M = 4 and 
N = 9, M = 6, respectively. The figures suggest that although (Y) depends on the 
precise distribution of the 5’s, this dependence is not very pronounced. Both cases result 
in (Y) which behaves like a sigmoid, namely it is stays very close to {Y){a = 0) = D 1 /D 2 
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Figure 2. The average yield (Y) as a function of cr for TV = 9, M = 6. The upper 
curve is for the Gaussian distribution (given by eq. (1^ 1 and the lower one for the 
bimodal distribution (given by eq. (12(1 b 


until it drops to zero over a relatively narrow range. This is indicative of the presence 
of a sharp crossover as a function of the noise level. 

The two hgures are not conclusive concerning whether for large a the ratio of the 
average yield obtained from the Gaussian distribution to that obtained from the bimodal 
distribution tends to 1, as predicted above (Figs. [T] and [2] only show that both tend to 
zero). Therefore, we present the ratio of the two for the case N = 6, M = 4 in Fig. [3] 
(the same hgure for the case iV = 9, M = 6 is quite similar and will not be presented 
here for brevity). 

It is of interest to obtain also the second moment of the yield. This calculation can 
also be performed analytically even for the Gaussian distribution but we will not present 
it in the following, because it is considerably more complicated and lengthy. In Figs. |4] 
- [5] we present the ratio of second moment to the square of the hrst moment (namely, 
{Y‘^) j fY)^) for the case N = M = 4 and the two distributions. Gonsidering Figs, 
m - 15] it becomes clear that the yield distribution behaves differently for three distinct 
regions of a. Indeed, that is the case as presented in Fig. |6] depicting yield distributions 
for three representative values of a. 

It follows from the above discussion that the average yield is a monotonically 
decreasing functions of a. It is actually also the case for its standard deviation 



demonstrated in Fig. [7] for iV = 6, M = 4, and can be 


proven more generally for the Bimodal and the Gaussian cases. 

At this point, a natural question that arises is whether the optimal yield is 
characterized by very small S’s? 
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Figure 3. The ratio between the average yield for the Gaussian distribution (given 
by eq. (l22ll 'l and that for the bimodal distribution (given by eq. (l20ll 'l as a function of 
cr for = 6, M = 4. 



Figure 4. The ratio between the second moment of the yield and the square of the 
mean / (Y) for the bimodal distribution as a function of cr for TV = 6, M = 4. 


To answer this question, consider the total energy (per period) 

N-M 

+ (23) 

m=0 

It is clear that the total energy is minimized when Bm = 0 for 0 < m < N — M and 
the minimal total energy per period is D 2 . We note further that the total energy can 
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Figure 5. The ratio between the second moment of the yield and the square of the 
mean / (F)^ for the Gaussian distribution as a function of cr for iV = 6, M = 4. 




Figure 6. The distribution of the Yield for three different representative values of a 
for iV = 6, M = 4: (a) a = 20, (b) cr = 65 and (c) a = 200, corresponding to the three 
regions identified in Fig. [5l As can be seen, the distributions are indeed very different 
with pronounced skewness in opposite directions. 


be expressed as E = y. Thus the (very simple) minimization process of E with respect 
to the choice of the B^s , can be expressed as 


1 dl I dY 

YdBZ~Y^dK 


(24) 


Consequently, on the set of 5’s for which the total energy per period is minimized we 
have 


dY 


dl 


dB^ ^ dB^ ’ 

and since Y is typically very small (in the superoscillatory regime), we conclude that 
at the point of minimal total energy, the yield is almost extremal (i.e., maximal) too. 
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Figure 7. The standard deviation of the yield cr(F) = y (Y^) — (Y)^ as a function of 
a for N = 6, M = A. As can be seen this is a monotonously decreasing function. 


because the absolute value of each of its derivatives is very small. Conversely, the 5’s 
at the optimal yield are very small, as expected from the a dependence of the average 
yield. Thus, optimization of the total energy under the constraints is typically almost 
equivalent to optimization of the yield itself. 

3. Noise in the Fourier coefficients of the signal 

So far, we have considered the behaviour of the yield in the space of interpolating signals, 
i.e. in the subspace where any choice of the i?’s yields a signal which automatically 
obey the constraints. Next we address the question of possible departure of the signal 
from the prescribed interpolation, in the vicinity of the yield optimized interpolating 
signal. Such departure will no doubt happen in any physical device that would need to 
utilize such signals, because any realization of a superoscillatory signal using a Fourier 
representation ([T]) would have access to the set of A’s with finite precision (due to the 
presence of noise for example). Such limited precision will obviously affect the yield, 
but of not less importance, it could also imply violation of the interpolation constraint, 
which may eventually destroy the superoscillatory nature of the signal. Recall that the 
yield is maximized with respect to the choice of the A^s under the constraints given by 
equations ([2])-(jl]). The result is a set of optimal Fourier coefficients |y4m|. To address 
the questions of violation of the constraints due to finite precision, assume that the 
actual Fourier coefficients describing the superoscillating function are given by 
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where are random errors and governed by a Ganssian distribntion with standard 
deviation 6. 

3.1. The interpolation sensitivity of the signal 

Consider the first qnestion, which is the sensitivity of the interpolation to random 
flnctnations aronnd the optimal solntion. Eqnations ©-(II]) imply that the change in 
the valne of the fnnction f{t) at the points tj = ja/{M — 1) for j = 0,M — 1, is given 
by 

N 

= (27) 

m=l 

where the absolnte valnes of all the Cjm’s are less or eqnal to 1 (by dehnition - see 
eq. ([2])). Conseqnently, 

\6fj\<xSy/N. (28) 

As long as the right hand side is (considerably) less than one we still have 

sign[f{tj)] = (-1)^, (29) 

so that the nnmber of oscillations in the interval [0, a] will not decrease and thus the 
required superoscillation is still preserved. 

To what extent do the fluctuations around the yield-optimized signal affect the 
required superoscillation? Put differently, what is the error level d that would start 
impairing the required superoscillations? From eqs. fl28|) - (l2^ it is clear that the required 
absolute accuracy is 

d<l/VN, (30) 

which does not seem to pose a serious problem. 

Note, however, that the problem might be that the error d can be reduced only below 
a value which is relative to the typical absolute value of the A’s. If that is the case it 
might pose a severe problem, because the optimal Fourier amplitudes are exponentially 
large in N (for superoscillating signals) and this is why out of the superoscillating interval 
the absolute value of / is exponentially large in N. In that case, the relative accuracy 
required in the Fourier coefficients is of the order oi 6 = d/ \ A\ oc where |A| 

denotes the typical size of the A^s and a > 1. Thus, if it is only relative accuracy that 
is possible to control, care should be taken that N is not too large in order to preserve 
the superoscillations. If absolute accuracy can be obtained, regardless of the magnitude 
of the optimal Fourier coefficients, it should be no problem to preserve superoscillations 
under very reasonable fluctuations in those quantities. We characterize the departure 
from the enforced interpolation, for a given set of errors in the Fourier coefficients by a 
single parameter, dV, which we call the interpolation sensitivity, and is dehned as 

.. M-l 

i=0 


( 31 ) 
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Figure 8. The mean interpolation sensitivity {5V) as a function of <5 for N = 6, M = A. 
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Figure 9. The mean interpolation sensitivity (SV) as a function of S tor N = 12, M = 
8 


Figures |8] - |9] show the interpolation sensitivity, which is the average of dV over the 
Gaussian distribution of the (5m’s, as a function of 5 for two examples of yield optimized 
superoscillations, 

Note that in spite of the fact that in the graphs above 5 attains values up to 10^, 
only for 5’s considerably smaller than one, we may expect the original frequency in the 
superoscillating interval to be unchanged. (The oscillation period, T, is taken as the 
average peak to peak distance and the frequency is a; = 27r/T.) 
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Figure 10. The normalized distribution, P{6V), for <5 = 0.01. The distribution is 
obtained from 10^ realizations of the set of Fourier coefficients. 


While the interpolation sensitivity (SV) is an important quantity that characterizes 
the deviations from the enforced interpolation, it is, obviously, not less important to 
obtain the distribution of the quantity SV. In the following we give this distribution for 
the case N = 12, M = 8 and for a specihc value of 5 = 0.01. It is clear that at least 
as far as the order of magnitude is considered, the interpolation sensitivity represents 
well the distribution of the values of SV. The distribution is obtained by considering 
10^ realizations of the set of the low frequency Fourier coefficients. 


3.2. Sensitivity of the Yield 


The next problem is that of the sensitivity of the superoscillation yield. The yield can 
be expressed in terms of the deviations as follows. The deviations from the optimal 
Fourier coefficients dehne a function giving the deviation from the optimal function, 
f{t), obeying the constraints, namely 


m = 


Si 


0 


N r 

— + ^ cos(mt). 

'' m=l ^ 


(32) 


The hrst order change in the yield can be written in terms of this function as 

/ f{t)S{t)dt f f{t)S{t)dt 

SY = 2'^ --• (33) 

/ Pit)dt J P{t)dt 

— TV —TV 


The fact that the yield has a linear correction when calculated at its maximum, should 
not be surprising, because it is a constrained maximum and the linear part appears as 
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Figure 11. The average deviation of the yield as a function of <5 for = 6, M = 4. 
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Figure 12. The average deviation of the yield as a function of 6 for N = 12, M = 8. 


a result of the small violations of the strict interpolation constraint. Had we considered 
deviations in the Fourier coefficients that respect the constraints, the lowest order 
correction would be quadratic in S{t). This implies that because of the departure from 
the interpolation constraints, the yield can increase under fluctuations. In Figs. ITTland 
[l2]we present the average of SY, as a function of 6. As we shall see, the average, which 
is not the average of the linear approximation of 6Y given above, is positive. This will 
be explained later when the scatter graph, giving the distribution of pairs {6V, 6Y) for 
given 6 will be discussed. 
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The relative deviation in the yield is given by 
/ f{t)d{t)dt J f{t)6{t)dt 


Y 


- 2 - 


/ p{t)dt 


I Pipdt 


(34) 


An upper bound on the relative error in the yield can be obtained by noting hrst that 
a t-independent bound on the mean square deviation function is easily obtained from 
equation (]3^ . 


71 

Then the error in the yield clearly obeys 

/ f{t)6{t)dt 


(35) 
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< 2 
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/ f{t)6{t)dt 


I PP)dt J P{t)dt 

—a —TT 

and by the Schwartz inequality and equation fl35l) above, 

aPN5 


(36) 
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TT 






a 

TT 

A 

1 

/ P{t)dt . 

/ P{t)dt 

— TT 


(37) 


The second term in the square brackets above can be neglected compared to the hrst 
and taking into account that within the superoscillating interval, P{t) is of order 1, we 
arrive at the conclusion that 




(38) 


where /? is a constant of order 1. 

Thus, if the condition for the preservation of superoscillation in the presence 
of additive errors in the optimal Fourier coefficients is obeyed, it ensures 

automatically a good upper bound on the error in the yield. In Figs. [13] - [T3| we 
show the average relative error in the yield, which we call the yield sensitivity. 

Since small changes in d result in large changes in the average of the deviation of 
the yield from its optimum (or equivalently to large deviations in the yield sensitivity) 
we may expect also the distribution of the yield to be very wide. The scale in the 
graphs describing those quantities suggest that if we would like to obtain a distribution, 
it would be reasonable to obtain the distribution of the logarithm of the yield rather 
than the distribution of the yield itself. Furthermore, after obtaining the distribution 
of 5V and the distribution of the logarithm of the yield, it would be interesting to ask 
if somehow the deviation from the constraints is related to the deviation in the yield. 
It is natural to assume that the interpolation constraint reduces the yield considerably. 
Therefore we may expect that the larger 5V is the larger is the yield. To check our hrst 
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Figure 13. The yield sensitivity (|^T|) /Y for N = 6, M = A. 


10 


> 10 

>- 

V 

10 


10 


10 


’ -2 -1 0 1 2 
10 10 10 10 10 

5 


Figure 14. The yield sensitivity (|5F|) /Y for N = 12, M = 8. 



hypothesis we present hrst the distribution of the logarithm of the yield (Fig. IThll . To 
check our second hypothesis, we present in Fig. [16] the scatter-plot of 6V and 6Y for a 
fixed 6. Namely, we chose a thousand sets {hm} from the Gaussian distribution with 
standard deviation 6. For each set we obtain 6V and 6Y and present the pair as a point 
in the 6V — 6Y plane. Both graphs (Figs. fThland fTHl) are presented for 6 = 0.01, iV = 12 
and M = 8. We see first that indeed, the distribution of the logarithm of the yield 
is narrow and resembles a Gaussian. We also see that increasing 6V tends to increase 
SY, with obvious tendency to have more positive dY's than negative ones. Therefore, 
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Figure 15. The normalized distribution of the logarithm of the yield, for 6 = 0.01. 
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Figure 16. The scatter-plot of {6V,6y) for 6 = 0.01, = 12 and M = 8. 

the initial insight that relaxing the constraints will result in higher yields seems to be 
correct. This is actually the reason why the average (SY) is always positive. 

3.3. Sensitivity of the number of superoscillations 

So far, we have characterized the departure from strict interpolation by the single 
parameter SV. It can occur, however, that while the departure from strict interpolation, 
as expressed by SV, is not that small, the signal still superoscillates with the prescribed 
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Figure 17. The average number of missing oscillations (An) as a function of of the 
error level <5 for = 15, M = 10. 


number of oscillation in the appropriate interval. The opposite is also possible. Namely, 
although 5V is relatively small, some oscillations may still be lost in the appropriate 
interval. With that in mind, we count the number of oscillations in the interval (—a, a), 
and denote the number of missing oscillations in this superoscillation interval by An (in 
all the examples we used a = 1). In Fig. [T7| we present the average of An as a function 
of 5 for the case A = 15, M = 10. As expected, the average depart from zero slowly, 
and stays fairly small even for (5’s that are not extremely small. The standard deviation 
of An as a function of 5 is also a useful tool for assessing the loss of oscillations as a 
function of the error level 5, and we therefore present it below in Fig. [181 

4. Summary and Conclusions 

To conclude, we have considered the effect of errors in the Fourier coefficients dehning the 
signal on the superoscillating yield and signal quality of interpolated superoscillations. 
Such errors are abundant in any physical system that may use superoscillations and 
therefore analyzing it is of high importance both theoretically and from the point of 
view of applications. 

We considered in detail two scenarios. In the hrst scenario, we analyzed deviations 
within the subspace of exactly interpolating signals from that signal that minimizes the 
signal energy per period. We have also shown that such a signal is very close to the signal 
of optimal yield within this subspace. The second scenario is where errors in the Fourier 
coefficients of the low frequency components can violate the interpolation condition and 
thus can reduce the signal quality in the superoscillation portion. Such errors affect 
also the yield. To quantify these effects we dehne the Interpolation-sensitivity 5V and 
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Figure 18. The standard deviation of the number of missing oscillations, as a function 
of S for N = 15, M = 10. 


Yield-sensitivity 6Y. Estimates on both sensitivities as a fnnction of the nnmber of the 
Fonrier components and the standard deviation describing the nncorrelated errors in 
each of the Fonrier coefficients have been obtained. We have also connted the nnmber 
of missing oscillations An dne to those errors. It is fonnd that an error level 6 in the 
Fonrier components {Am} which is not too small still ensnres snperoscillation and can 
resnlt in a yield which is larger than the optimal one nnder the interpolation constraint. 
Namely, a slight relaxation of the strict constraint may still keep the reqnired nnmber 
of snperoscillations yet resnlt in an increased yield which can be a real advantage in 
application. 

An important bottom line is that althongh generating optimized snperoscillations 
may be a challenging task, as it typically reqnires very high machine precision, once 
attained they are no longer very sensitive to error or noise, and can thus be realized 
in physical systems where noise is abundant. Another important conclusion is that 
compromising the shape of the function in the superoscillatory region (such as the 
regularity of the oscillations and the location of their maxima/minima) can result in an 
improved superoscillatory yield, which could be advantageous. 
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